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f^ The effect of scattering with non-zero impact parameters between consituents in relativistic heavy ion colUsions is investigated. 

T-H : In solving the relativistic Boltzmann equation, the characteristic range of the collision kernel is varied from approximately one 

^^ fm to zero while leaving the mean-free path unchanged. Modifying this range is shown to significantly afi^ect spectra and flow 

(-H , observables. The flnite range is shown to provide effective viscosities, shear, bulk viscosity and heat conductivity, with the 

■*-^ ■ viscous coefficients being proportional to the square of the interaction range. 



I. INTRODUCTION 



A principal goal of relativistic heavy ion collisions is to experimentally discern bulk properties of the excited vacuum. 
Si^ ■ To accomplish this aim it is imperative that one understands the implications of the finite size and lifetime of the 
?H ■ global reaction. Two microscopic length scales govern the importance of finite-size effects, the mean free path and the 
range of interaction. Microscopic models, e.g., those based on the Boltzmann equation, easily incorporate the effects 
of a finite mean free path. Such effects can be linked to viscous terms in analogous hydrodynamic descriptions. 

In intermediate-energy collisions, where excitation energies are tens of MeV per nucleon, the role of the finite range 
to the strong interaction has been studied in its relation to the surface energy of nuclear matter. In such Boltzmann 
descriptions the binding energy of nuclear matter is introduced via the mean field, with the coarseness of the mean field 
mesh being adjusted so that the effect of the interaction range is effectively tuned to reproduce the surface energy 
of nuclear matter. In non-relativistic molecular dynamics, the effects of hard-sphere interactions have also been 
investigated. In this case the size of the spheres represents a length scale which can strongly affect bulk properties 
of the matter at high density. However, in the context of a Boltzmann description, where n-body correlations are 
explicitly neglected, the effects of a finite range inherent to the scattering kernel have not been analyzed for their 
impact on final-state observables. 

It is the goal of this study to ascertain the importance of this second length scale in high-energy collisions. By 
varying the interaction range in the scattering kernel, while leaving the mean free path unchanged, we study the 
manisfestations of a non-zero interaction range. We demonstrate that a finite interaction range contributes to viscous 
terms in a manner similar to the finite mean free path, but with different dependencies with respect to density and 
temperature. We find that the finite range of the interaction affects the evolution of heavy-ion reactions, and alters 
final-state observables, especially the elliptic flow. 

In the next section we present a formal review of the Boltzmann equation and show how viscosities arise from 



the interaction range. In section HI we show how spectra and flow observables are sensitive to the finite range. We 



discuss algorithmic sensitivities which cannot be ignored when the interaction range is non-zero, and compare results 



compiled from four similar numerical implementations of the Boltzmann equation to illustrate this sensitivity. We 
show how viscous heating can be understood from analyzing the collision kernel in Sec. KM. In particular we present 
a comparison of heating derived from analysis of the collision kernel with the heating observed in a simple simulation. 
A discussion of the importance of non-local interactions and the associated causality problems is provided in the 
summary. 

II. CONNECTING VISCOSITIES TO FINITE-RANGE INTERACTIONS 
A. The Role of the ColHsion Kernel in Boltzmann Descriptions 

The Boltzmann equation can be expressed, 

|/,p,M) + V, . V/(p..O +F,.0 . V,/(p,.,,) = /.=,.=,VV.V.,' ,1) 

• {/(q, r, i)/(q', r', t')K{r -v\t-t'- q, q'; p, p') - /(p, r, t)/(p', r', t')K{v -r',t- t'; p, p'; q, q')} , (2) 

where / is the phase space density and F is the force dp/dt felt by a particle at position r. The collision kernel 
/C(r — r', t — t'; q, q'; p, p') describes the differential probability for scattering a pair of particles separated in space- 
time by X — x' with initial momenta q and q' into final states with momenta p and p'. The range of the collision 
kernel in coordinate space is the subject of this paper. 

Integrating over the collision kernel should yield the cross section, 

/I /-7 

d^r'dt'Kir - r', t - i'; q, q'; p, p') = _-_^^z;,,i<S3(p + p' - q - q'), (3) 

where prci is the relative momentum of the outgoing pair (p — p')/2. By inspection of Eq. (||), one can see that 
the coordinate-space dependence of JC appears rather arbitrary as long as it integrates to the free-space cross section. 
Indeed, results at low density, where particles interact only pairwise, are unaffected by the form of /C as long as the 
range is much less than the mean free path and much less than the characteristic dimensions of the reaction volume. 
Any non-zero extent of the collision kernel leads to problems with super-luminar transport. However, these problems 
are easily defeated by restricting the kernel to being local, i.e., 

K{r -r',t- t'; q, q'; p, p') = j^S^r - r')S{t - t') ■ -^v,,,d^p + p' - q - q')- (4) 
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The Boltzmann equation can now be written in a manifestly covariant form. 
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Here, Up is the four- velocity of a particle with momentum p, F^ is the force dp^ /dr, and Prci is the relative momentum 
of the outgoing particles in the center of mass. 

B. Effective Viscosities from Finite-Range Interactions 

In this section we describe how interaction over a finite range contributes to the shear viscosity, 77, the bulk viscosity, 
C, and the thermal conductivity, x- We relate the range of the interaction to all three coefficients. In order to make this 
connection, we consider two particles which scatter from one another, separated by a distance r = r2 — ri. Combining 
this finite separation with the velocity gradient, one sees that the first particle interacts with particles which have a 
higher average energy. By evaluating the rate at which energy is transferred to the first particle from colliding with 
more energetic particles, we find an expression for the rate at which heat is deposited to the region defined by ri. By 
comparing to analogous expressions from hydrodynamics, we derive expressions for all three viscous terms in terms 
of the interaction range, r = |ri — r2|, the density n, and the collision rate T. 



Choosing a reference frame such that the velocity of bulk matter at the location of the first particle is zero, the 
collective velocity at r2 is 



Aij Tj , Aij 



dr-i 



(6) 



For an elastic collision where two particles of identical mass simultaneously change their momenta, the radial 
components of the momenta must be interchanged by the collision if energy, linear momentum and angular momentum 
are to be conserved. Physically, this corresponds to the scattering from the interior or exterior surface of a hard sphere. 
The average energy change of the first particle is then 



{5E^)^{E2,r)~{E^,r) = ^{{^-ff). 



(7) 



The mass term m is not to be taken literally as the mass of the particles, since the averaging may include factors of 
the velocity to account for the flux or it may have a complicated form to accommodate a desired differential cross 
section. For relativistic motion, the mass might incorporate the lateral motion of the particles. Writing (SEi) in 
terms of A, 



{SEi 



m 
2^ 
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(8) 



One can perform the average over the directions of r using the identity, 



AjAkiinrjrkri) = AijAki— {SijSu + SikSji + SuSjk) ■ 



One can then express {5Ei) in terms of A and r. 
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The rate at which the entropy increases due to these interactions is given by the density multiplied by the rate at 
which collisions deposit energy non-locally, 
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Here, F is the collision rate experienced by a single particle and n is the density. 

It is notable that only the symmetric part of A contributes to {6Ei). This owes itself to conservation of angular 
momentum which forbids rotational motion from being transferred between particles. In fact, if one had derived an 
expression for (SEi) using w^ instead of v'^, the resulting expression would have included the odd parts of A which 
would correspond to rotational motion, V x v. These terms would have no hydrodynamical analog as they would 
have reflected a violation of angular momentum conservation. 

We now provide analogous expressions for Eq. (11) in the language of hydrodynamics. The expression for entropy 
production ffl in terms of velocity gradients is 



^^«^^E(£ + fe-|*«V^v)4,V^v,».i,VT,^ 



(12) 



One can extract the coefficients by comparing Eq. (^Tj) to the first two terms in Eq. (p^). 



mnTr^ 

C = ^. (14) 

Applying similar reasoning, one can also derive an expression for the thermal conductivity. First, we relate the 
temperature gradient to the energy flow. Again, we consider particles separated by r. If collisions occur between two 
particles at locations with different temperatures, the average energy exchanged is 

SEi = ^CrT ■ VT. (15) 

Here Cr represents the change in radial kinetic energy per particle per change in temperature, 

a = ^Er. (16) 

In the non-relativistic limit, Cr = 1/2. 

Since the exchange corresponds to moving an energy a finite distance over an effective time given by the collision 
rate, one can define the average momentum density in terms of the energy flow. 

nTCr, .dT 
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An extra factor of 1/2 was added to correct for double counting the collisions. 
One can relate the energy flow to the entropy production, 

dS f ,. Ide 
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d'-'^^i^Tr, (21) 
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where the continuity equation has been applied. One can now compare to the last term in Eq. dl2) to obtain the 
thermal conductivity x- 

(22) 
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The forms for the three coefficients, ry, C and x-, fundamentally differ from the forms that result from considering 
a non-zero mean free path £. The viscous coefficients that result from finite i are independent of density and scale 
inversely with the cross section |1|. The coefficients arising from a non-zero interaction range r scale as the square 
of the density and are proportional to the cross section. Thus, non-local interactions provide viscosities that are 
important in wholly different conditions than those where the finite mean free path plays an important role. For a 
rapidly expanding system, finite-range interactions play an important role when the interaction range multiplied by 
the velocity gradient provides velocities of similar magnitudes to local thermal velocities. Such conditions exist in the 
first one fm/c of highly relativistic hadronic collisions. 

The expressions here derive from a very specific picture, non-relativistic particles moving on straight-line trajectories 
punctuated by sharp collisions when the separation equals r. However, all interaction at a finite distance should result 
in viscous behavior. Relating the distance r to the cross section might involve a detailed microscopic evaluation of 
the collision kernel. This is especially true for relativistic motion. Despite the complications, one indeed expects the 
distance r^ to be of similar magnitude to collision cross sections. It would be interesting to discern whether the ratio 
between the shear and bulk viscosities varies for different scattering models. 



III. SPECTRA AND FLOW 
A. Sampling Factors and Scattering Algorithms 

Numerical realization of the Boltzmann equation is usually accomplished by performing a classical simulation of the 
event and applying a large over-sampling factor A PJ3]. For instance, if one is sampling a system with 100 particles 
the Boltzmann simulation might simulate the trajectories of 3200 particles which would correspond to A = 32. The 
scattering kernel must be diminished by a factor of 1/A to ensure that the scattering rate and mean free path is 
independent of A. As the Boltzmann equation is built on the hypothesis that the evolution is completely described by 
the one-body distribution, the Boltzmann limit is realized for large A as the correlations between particles becomes 
negligible. 

The most common way in which the scattering rates are scaled by 1/A is through reducing the cross section. 
Usually, particles that pass within a distance a/ct/tt of one another are scattered. Scaling the cross section by a 
factor of 1/A leaves the collision rate and mean free path unchanged, and in the limit of large A leads to a purely 
local scattering kernel which solves numerous problems associated with non-causal transmission of information and 
invariance to Loretz boosts. 

Although one can solve the covariance problems by choosing a large A, one must then understand whether ignoring 
the finite interaction range has affected the results. More realistically, particles do interact at a finite range through 
mutually coupling to classical fields or through the quantum exchange of particles. A variety of ideas have been 
discussed in the literature for including such interactions in ways which preserve covariance and are consistent with 
the uncertainty principle H. However, due to an assortment of technical challenges, all such approaches remain in the 
development stage. Since most comparisons with experiment continue to be performed with Boltzmann approaches 
and since understanding the effect of finite interaction will be more clear with a simple model, we begin our study of 
finite-range effects by analyzing classical Boltzmann treatments. 

For the investigations described in this subsection, we will always associate large sampling factors with a reduction 
in the effective range of the interaction. However, that is not necessarily the case. Instead of reducing the spatial 
extent of the cross section, one could reduce the probability of scattering by a factor of 1/A, i.e. not all particles 
passing within a distance \/oJtt of one another would scatter. For high sampling rates collision-finding meshes with 
granularities determined by the interaction range are necessary for efficient numerical algorithms. By reducing the 
cross section, instead of reducing the scattering probability, one is able to reduce the mesh size which results in 
significantly faster execution. 

As mentioned above, particles typically scatter in Boltzmann algorithms when they pass within yapK of one 
another, at a time chosen to correspond to the point of closest approach. Aside from choosing the range of the 
collision kernel, many other aspects of such algorithms are arbitrary. For instance, one might expect that particles 
should scatter with an impact-parameter-dependent probability that has a more complicated form than a simple step 
function. Even in the non-relativistic limit, scattering at the point of closest approach violates conservation of angular 
momentum. In some algorithms the reaction plane is preserved in each scattering which has been shown to affect 
resulting flows in the cascade limit (A = 1) J5l. Finally, scatterings need not occur instantly. By delaying the particles 
a certain time before they emerge with their asymptotic momenta, one can effectively alter the pressure. One must 
coordinate the time delays with in- medium modifications and knowledge of the phase shifts in order to be ergodically 
consistent with thermodynamic properties at the level of the second virial coefficient [|6|] . 

In the cascade limit, all scattering algorithms carry a degree of arbitrariness due to issues involving Lorentz invari- 
ance Q:^]- One problem deals with the ordering of collisions in time. Each collision can be assigned a point in space 
time. Since many points would have space-like separations, time ordering is frame dependent. Finally, since each 
collision involves changing trajectories over a finite separation (when A ^ oo) the final outcome can be affected by 
the time-ordering prescription. 

A second aspect of arbitrariness deriving from relativistic considerations derives with what is meant by "point of 
closest approach" . Since altering the two trajectories simultaneously is inherently frame dependent one must choose 
a prescription for the time at which each particle is to change course due to the collision. To illustrate this we 
consider two particles with momenta pi and p2 at space-time coordinates xi and X2- One can solve for the points 
on their trajectories, x\ and x'2^ where they would reach the point of closest approach as measured by an observer in 
the two-particle center-of-mass, by solving for the point where the relative separation becomes perpendicular to the 
relative momentum, 

^'.M^^M+pM^M^_^ (23) 

Pi ■ r" 



o^T = <+p' /'f~'''^ (24) 

0^{x[~x',)-{q~P^). (25) 

Here P and q are the total and relative momenta respectively, and P ■ X/^/s is the time at which the collision occurs 
as measured by an observer in the center-of-mass frame. One can solve Eq.s ( p3[j25| ) for P ■ X , then inserting into 
Eq.s ( p3| ) and (|23) one can solve for x[ and Xj. 

An algorithmic choice that would retain covariance at the two-body level would be to alter the two momenta at x'l 
and x'2 even though the times would not be simultaneous in other frames. However, in a dense medium particles may 
be involved in several collisions which may result in a variety of events being scheduled between the times t[ and tj- 
In fact, a particle might be assigned a collision time that precedes its creation. Additionally, Boltzmann algorithms 
often incorporate collision-finding or mean-field meshes. Prescriptions where the two scattered particles have their 
trajectories altered at different times become problematic when coordinating the evolution with the mesh. The most 
common algorithmic choice involves choosing a fixed frame to define collision times, then altering the trajectories 
at the point of closest approach as measured in that frame. This prescription becomes questionable when collective 
velocities become relativistic. Such is the case at RHIC where particle emission is spread over 10 units of rapidity. 
This problem can be somewhat alleviated by ordering the collisions according to r = (t^ — z^)^/^. This choice seems 
reasonable for ultra-relativistic collisions where conditions are determined largely by t, but seems unreasonable for 
lower energy collisions where longitudinal boost invariancc is not realized. 

B. Comparing Results from Four Algorithms 

To illustrate the sensitivity of results to algorithmic choices, we compare results from four numerical approximations 
to the Boltzmann equation which are realized in four separate models, each arising from a different author. All the 
algorithms discussed here scatter particles instantaneously at the point of closest approach, and all the models assume 
a simple s— wave form to the cross section. Each algorithm was executed from the same initial conditions, and each 
was executed both in the cascade limit (A = 1) and with a high sampling factor, A = 16. Due to the arbitrariness 
inherent to algorithms with finite-range interactions it should not be surprising for results to vary between codes with 
A = 1, but it is expected that results from the four approaches should be indistinguishable in the large A limit. 

The four models differ in their definition of "simultaneous" , in their time-ordering prescriptions and in whether they 
preserve the scattering plane in two-particle scattering. The four algorithms compared here represent the following 
choices: 

1. ZPC: A parton cascade code authored by B. Zhang ||]. Particles scatter at the point of closest approach, 
simultaneous in a fixed reference frame. This time is determined by first finding the space-time points x'l 
and x'2, where the two particles would scatter should they alter their trajectories simultaneously in the two- 
particle center-of-mass frame. In a designated laboratory frame, a time is determined by averaging t'l and 
^2. The trajectories of both particles are simultaneously altered at this average time. The reaction plane is 
preserved in each two-body scattering, and rescattering between the same partners is allowed after one particle 
has rescattered. 

2. GROMIT-i: A generic scattering engine developed for the RHIC Transport Theory Collaboration (RTTC) [|o). 
This should be identical to ZPC. 

3. MPC(0.4.0): A Boltzmann description authored by D. Molnar ||ll|. The principal difference with the codes above 
is that the reaction plane is not preserved in the two-body scattering. Another difference is that rescattering 
between collision partners is not allowed until both particles have rescattered. It should be noted that A/IPC 
has switches which make it possible to reproduce the choices described in ZPC and GROMIT-t. 

4. GROMIT-t: Another version of cascadc/Boltzmann engine developed by the RTTC collaboration with collisions 
ordered by r to be appropriate for ultra-relativistic collisions |10| . The space-time points at which the two 
particles would scatter had the particles scattered simultaneously in the two-particle rest frame, x'l and x'2, are 
used to generate two proper times, r{ and T2. The average t' — (t{ + T2)/2 is then used for ordering. Both 
trajectories are altered when their proper time equals r'. The reaction plane is preserved in the scattering and 
scattering between the same pair of particles is allowed after one of the particles has rescattered. 



Calculations were performed for all four models with identical initial conditions based on a thermal Bjorken geometry 
L2[. In Bjorken coordinates, the time is replaced by a proper time, 



T 



= ^t^ - z^, (26) 

and the position along the beam axis, z, is replaced by the effective coordinate, 

,4.o.(i±f). (... 

The coordinate 77 also gives the rapidity of an observer moving at constant velocity from the origin at t = to the 
space-time point {t, z). 

Thermal emission from a Bjorken geometry would be realized by choosing thermal sources with positions equally 
distributed in 77, with source rapidities j/s = rj, and with emission taking place at a fixed proper time, r. Such emission 
would appear invariant to boosts along the beam axis. In the calculations described here, the sources are confined to 
a finite range in 77, 

- 2 < ?7 < 2. (28) 

The transverse spatial coordinates were randomly placed in a cylinder of radius 5 fm. The simulation involved 2400 
pions and 240 protons which were assigned momenta according to a thermal distribution with a temperature T = 180 
MeV for the protons and T = 165 MeV for the pions, to roughly acheive consistency with spectra resulting from heavy 
ion collisions. The initial particle densities were 7.64 pions and 0.764 protons per fm^^. To be as simple as possible, 
constant cross sections were imposed independent of the species involved. 

The resulting spectra at mid-rapidity as calculated with GROMIT-t are shown in Fig. |l| for four different cross 
sections, 0, 10, 20 and 40 mb. These calculations were performed with a sampling factor A — 32. For larger cross 
sections, the pion spectra are cooler while the proton spectra are somewhat hotter. The difference would be more 
pronounced had the pion and proton distributions been initialized with identical temperatures. 

Figure displays spectra from GROMIT-r with a =40 mb for four different sampling factors, A = 1, 2, 8 and 32. 
The sensitivity to A is remarkable. For small A it appears that both the pion and proton spectra become hotter. 
Results appear to converge for A ^ cx) as the A = 8 and A = 32 results are barely distinguishable. In all the 
models investigated here, the final-state transverse energies were higher for small sampling factors. We attribute this 
reduction in cooling to viscous effects arising from finite-range interactions. 

In the upper panel of Fig. g the resulting spectra from the four algorithms are displayed in the cascade limit. The 
slopes vary by approximately 20 MeV. The same results are displayed in the lower panel of Fig. but with a high 
sampling factor, A = 16. The various algorithms are then in agreement. Since the hadrons are initially separated by 
approximately 1/3 fm, while the interaction range is 1.1 fm, it is not surprising that the A = 1 result is sensitive to 
the scattering algorithm. 

For each algorithm the spectra fall significantly more steeply when produced with a higher sampling factor A. This 
is especially true for the proton spectra, which suggests that radial fiow might be stronger in the cascade limit. Since 
the average transverse energy for pions is also larger in the cascade limit, one can infer that longitudinal cooling is 
significantly suppressed through particles interacting at a finite range. 

C. Elliptic Flow 

As a second exploration of the sensitivities to the sampling factor, we consider elliptic fiow as observed from the 
GROMIT-r algorithm. In this case cyclic boundary conditions in rj were employed to simulate a truly boost-invariant 
system. The densities of protons and pions were chosen identically as above, 7^,r = 7.64 fm~'^, Up = 0.764 fm~^ at 
T = 1.0 fm/c. The temperatures for the pions and protons were again chosen to be 165 MeV and 180 MeV respectively. 
Aside from the cyclic boundary conditions in 77 the only difference with the initial conditions used in the previous 
subsection is that the particles were confined to an ellipse rather than a circle in the transverse direction. The two 
radii of the ellipse were 5 fm and 2.5 fm. 

The elliptical shape resulted in elliptic fiow which is parameterized by the observable V2- 

7;2 = (cos(0-0o)), (29) 

where (f>Q specifies the direction of the reaction plane which contains the short axis of the ellipse. Elliptic flow has 
been proposed as a means for determining the equation of state of the matter at early times H-EO| . 



Resulting flows are shown in Fig. ^ and Fig. |5| for two cross sections, 10 mb and 40 mb, and for four sampling 
factors. In Fig. H the cross section was left constant as A was varied. Instead, only a fraction, 1/A, of the possible 
scatterings were executed to account for the over-sampling. The lack of a sensitivity to A suggests that the Boltzmann 
limit is reached even for A = 1, i.e. n-body correlations do not signficantly affect the flow. The same quantities were 
calculated in Fig. ^, except that the cross sections were scaled as 1/A to account for the over-sampling. The results 
are then highly sensitive to A which suggests that the finite interaction range that accompanies small A leads to a 
reduced elliptic flow. 

The reduction in longitudinal work and the damping of elliptic flow due to the incorporation of a finite-range 
interaction into a Boltzmann treatment suggests that viscous effects have effectively been introduced into the reaction. 

IV. ANALYZING THE COLLISION KERNEL 



In Sec. II B a finite range of interaction, characterized by a length scale, r, was shown to generate viscous behavior. 
This length scale is determined by the cross section, r^ r^ ct/tt, but the constant of proportionality is not trivially 
determined and can depend on seemingly arbitrary aspects of scattering algorithms. In the following subsection, we 
illustrate how viscous heating can also be directly related to the collision kernel by considering a simple example of 
one-dimensional expansion with a Bjorken space-time geometry. In the subsequent subsection we compare predictions 
based on the form of the scattering kernel with results from simulations based on the same kernel. 

A. Viscous Heating in a Bjorken Expansion 

As a simple example we again consider a one-dimensionally boost-invariant system of infinite extent in all dimen- 
sions. We keep the number of particles fixed so that the density scales inversely with r, 

„(t) ^ -dN/{dAdr]). (30) 

r 

where the number of particles per area per rapidity interval, dN/{dAdri), is fixed. To simplify matters we consider 
massless particles. 

Here we wish to calculate the rate at which a particle at r/ = has collisions with other particles, dNc/dr, and the 
rate at which it gains or loses energy from such collisions, d{Ei)/dT. We assume the particles are in local thermal 
equilibrium, and analyze the collision kernel to determine the rates. 

Refering to the colliding partners with the subscript '2', the collision and heating rates per particle are 

dNc dN f ^ ^ . dN 



■S{r^-r)e[^-b'), 



d(Ei) dN f ^ dN 

= 271— — / dr,2r2dr2d^P2:7^ (32) 



dr dA drj J d^p2 

Here, SEi is the average change in energy of the particle '1' due to the collision and r2 and 772 describe the position 
of the second particle. The impact parameter for the two-particle collision is b and Tc is the collision time. 

The strategy employed here is to calculate the collision time, Tc- The first particle is at a position xi = (r, 0, 0, 0) 
with four momentum pi. The second particle is at a position X2 — (TC0sh772, '^2: Oi''"sinh772), with a four- momentum 
P2- Once Tc is understood in terms of r2, 772, Pi and p2, one can replace the delta function in the expression above, 

S{Tc-T)-.S{rc~r2)j-^, (33) 

where r^ is the position required to make the collision occur at r. By substituting the delta function with Vc for the 
delta function with Tc, the integrals in Eq. (pl^) and Eq. ( p2| ) can be simplified and solved numerically. 

The first step one must perform is to find Tc in terms of r. The prescription for finding Tc is somewhat arbitrary 
due to covariance issues. For our purposes, we will first define a time in the rest frame of the two particles where 



the particles run abreast of one another, tc- With Xc — (tc, 0,0,0) in the center-of-mass frame and P = pi + P2, the 
conision time is determined by the condition, 

(pi-p2)-(x'i-x^)==0, (34) 

Xi=xiH pi, (35) 



One can solve for P ■ Xn 
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One can now solve for x[ and x'2. 



X2=X2^ 5 P2- (36) 

P ■P2 

P ■ Xc^ Pi- XI+P2- X2. (37) 

' , P2-{X2- Xi) 

Xn^xi+pi (38) 

Pi ■P2 

, , Pi ■ {xi -X2) , . 

X, = X2 + P2 (39) 

Pi ■P2 

If one calculated the times at which the collision occurred, t'l and t2, they would only be simultaneous in the two- 
particle center-of-mass frame. Similarly, if one calculated t{ = {t'^ — z'-^Y^^ and T2 = {t'2 — z'^Y^"^ ^ one would find 
that the two Bjorken proper times would not be simultaneous. One must arbitrarily choose a prescription to choose 
Tc in terms of x'-^ and x^- For our purposes, we choose the following prescription, 

rl ^ ^±^ (40) 

Our only motivation in averaging the r^s rather than the ts is that the algebra simplified. Using this choice, 

rl . ,^ + C,r-^2Cir'^C2 ^^^ 

C2 2 1 2 2 
== 7lPl_L + 72P2_L 

Cl == 71^1 + <5l7lP?_L + 72^2 + hl2p\\_ 

C2 = ^\p\y_ + 2(5iai + (5^P2_L + 2<52a2 
The coefficients are defined 

a\=Ext-pY^y_z ,a2=E2t + P2,7.^ (42) 

(5i = 2zp2.z/{pi ■P2) ,(^2 = -2zpi^:z/{pi ■P2) 

71 = '^P2.x/{Pl ■P2) ,72 = -2pi,a;/(pi ■P2) 

Here t, r and z are the positions describing the first particle in a frame where the particles are centered about 77 = 

and r — Q. 

V2 T2 

rl ~ — , t~ T cosh 77', z — — Tsinhyy', r' ~ — . (43) 

In terms of these new variables, one can rewrite the rates above, 

dNr _,, dN f , , , , , ,^ dN 

16TT-——dr]'r'dr'd^P2-j:^— (44) 



dr dA drj J d^p2 

•<5(r-r,)e(- 
Vtt 

= 167r-— -— / drj' r' dr' d^p2 



dA drj J d^p2 

-^<"'--'ft;^«(^'1- 

d(Ei) dN f , , , , , ,^ dN 

= l^TT— -— / dr]' r' dr' d^p2-:^— (45) 



dr dA drj J d^p2 

■5{r'-rc) ^ f^ 9(^-6^) ^i^i. 
Lor + Gi Vtt / 



Core + Ci 


/ 


Vtt 


^b' 


re 


SEi 


Cor, + Ci 


Vtt 


< 



When solving for r,, there are two solutions to Eq. ([4l|). If one uses solutions for both positive and negative r,, the 
factor IGtt should be reduced to Stt. 

One can express these results as averages over pi and p2 



If scattering angles are chosen with equal probability forwards and backwards, the average change in the two energies 
is 

SE, = ^^. (48) 

The non-local aspect of the collision kernel should contribute to heating the particles in their local frame. From 
physical arguments, one expects the non-local contribution to heating to scale with temperature, density, time and 
cross section in a simple manner. 

dEi_ ^ _3_d(E^ ^^g^ 

dr An dr 

dN a'^T 



= /? 



dAd-q 



where Et is the mean transverse energy per particle and /3 is a dimensionless constant. 

The simple scaling derives from the hydrodynamically motivated Eq.s (12) and (|l3|). One expects the collision 



rate to scale proportional to the density, which requires the factor dN / {dAdrf} and one inverse power of r, and 
the cross-section. The squared- velocity gradient suggests an extra factor of r~^, and the range of the interaction 
requires an extra factor of a. The constant j3 is determined by the form of the differential cross section, e.g., s-wave 
scattering would result in a higher f3 than a highly forward-peaked form. Since one power of a comes from the range 
of interaction, (3 should scale inversely with the sampling factor A. Finally, the effective mass should be proportional 
to the temperature. 

The heating rate due to non-local interactions is calculated from Eq. ( |47| ) by numerically analyzing the collision 
kernel and is displayed in Fig. ^ after being scaled by the temperature, time and cross section. Had the simple scaling 
argument been correct the ratio would have been a constant /3. However, due to higher order corrections in l/r, 
the ratio varies as a function of t. The ratio approaches a constant for large r, r^ >> a. The scaling fails when r 
becomes less than the range of the interaction, which for this example is ~ 1.0 fm due to the 10 mb cross section. The 
departure of the ratio from a straight line illustrates the limitation of simple hydrodynamic arguments to describe 
the viscous heating from non-local interactions. 

B. Comparing to Numerical Results 

Finally, we present numerical results involving a boost-invariant Boltzmann description (GROMIT-r) which was 
executed with cyclic boundary conditions, both in rj and in the transverse coordinates. We compare the viscous 
heating rate observed in the numerical calculation with the viscous heating rate expected from the scaling arguments 
expressed in Eq. ([49|), where the coefficient (3 was determined from analyzing the collision kernel from the last section. 
The temperature was set to 400 MeV at a time r = 0.1 fm/c, and the cross section was chosen to be 10 mb. A simple 
s-wave form was used for the angular dependence of the cross section. 

The resulting mean transverse energy is displayed in Fig. H as a function of r. The initial heating derives from 
the non-local nature of the interactions. Longitudinal cooling ultimately dominates the behavior as the non-local 
contribution to the heating falls roughly at r^^. 

The dotted line in Fig. |^ describes the evolution of the transverse energy in the limit of ideal (non-viscous) 
hydrodynamics. In that limit, the stress-energy tensor has a simple form, 
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T"'' = ew"M^ + P(u"'u'3-g"'3). (50) 

In the Bjorken limit, dv/dz = \/t and the evolution of the energy density becomes 

For the massless case, P — e/3, which gives the result, 

e{r) = e{To) (^- j (52) 

The velocity gradient generates a shear which contributes an additional term to the stress-energy tensor. 



T,-^,(-.--»"«T(||^ + |»J<a-.)^ (53, 

In a simulation, one can determine ryjvs by evaluating the stress-energy tensor. 

_T^^^2^^. (54) 



T^^'+Tyy ^,, ^VNS 



The components of the stress-energy tensor are extracted by analyzing the momenta of particles as measured in the 
local rest frame. 

The Navier-Stokes evolution of the energy density is governed by the equation, 

OT T 6 T^ 

For massless particles interacting with a constant cross section, dimensional arguments force the viscosity to rise 
linearly with r since the mean free path should grow with r due to the density falling as 1/r. 

TINS = Cns^t, (57) 

where C^s is a constant, determined only by the form of the cross section. The energy density then follows the Navier 
Stokes form pl|. This form should be valid unless the viscous contribution to the stress energy tensor approaches 
the same order as the equilibrated pressure p^ ]. 

/^\ -(4/3)(l-Civs) 

e(r) = e{r) ^-j . (58) 

The Navier Stokes result is represented by a dashed line in Fig. 0. The value of Cats was determined by evaluating 
the asymmetry of the stress-energy tensor in the simulation at large times. Running simulations with a large sampling 
ratios generated results in excellent agreement with the Navier Stokes result. 

The inclusion of non-local effects is responsible for the discrepancy between the simulation results in Fig. and 
the Navier-Stokes results. Eq. ( p^ ) suggests that the non-local correction to the viscosity should scale proportional 

tOT"^ 

Vni = Cnie/r, (59) 

where C„; is independent of r and scales with a^ as explained in the previous section. One can now determine the 
evolution of the energy density by solving the hydrodynamical equations of motion, 

9 4 e 4 e 

g^e^--il-C^s)- + ^C^i^ (60) 

(61) 
This equation has a simple solution. 

11 



4/3(1-Cns) ( 2(J 



This form is shown with the sohd hnes in Fig. 0. The constant C„i was determined by the asymptotic hmit of Fig. 
y. The effect of non-local interactions is somewhat overestimated at small times by Eq. (62) as would be expected by 
considering Fig. || which shows that that growth of the viscous heating at small r is somewhat slower than the naive 
expectation that it scales as r^^. 

The effect of non-local interactions is lessened for scatterings that are more forward peaked. Fig. g illustrates the 
behavior of Et as a function of r for a screened Rutherford scattering, 

da 97r al 

(63) 



dt 2 {t-^i'^f 
97ra? 



(64) 



2Ai2 ' 
t={pi-p^f. (65) 

Here, as = y2/9, and the screening mass ^ is chosen to provide a cross section of 10 mb. As compared to the 
s-wave scattering result in Fig. M, the effect of non-local interactions is diminished while the Navier-Stokes viscosity 
is increased. This is expected since the mean free path is effectively increased while the energy transfer inherent to 
collisions is decreased. 

V. DISCUSSION AND SUMMARY 

The effects of short-range, but non-local, interactions provide an unwelcome complication to the analysis and 
interpretation of experimental results from relativistic heavy ion collisions. Flow observables are most significantly 
affected, especially elliptic flow. Thus, the additional length scale, characteristic of hard collisions, interferes with our 
ability to extract bulk matter properties due to the finite volume in space-time subtended by heavy ion collisions. 

By increasing the sampling ratio in Boltzmann treatments while simultaneously reducing the cross section, one can 
effectively eliminate these effects in models. Aside from simplifying the analysis, large sampling ratios are attractive 
as they eliminate sensitivities to a variety of arbitrary choices inherent to simulations and solve a variety of problems 
related to acausal propagation. However, in nature particles indeed interact over a finite range, either by exchange of 
off-shell particles or through a mutual interaction through classical fields. Therefore, it is important to understand 
the degree to which these effects are physical as opposed to representing numeric artifacts. 

Determining the correct range for non-local interactions in the first fm/c of a collision might require correctly taking 
into account both relativity and quantum effects. For instance, one might determine the extent by considering the 
transfer of momentum and energy through exchange of off-shell particles in the context of a Wigner decomposition 
of the propagator M. This object might then be analyzed in a matter analogous to the method in which the classical 



scattering kernel was studied in Sec. IV. However, this is beyond the scope of the current paper. 



On the optimistic side, non-local effects can be interpreted in terms of viscous parameters as shown in Sec. HB. 
This allows the incorporation of non-local effects into hydrodynamic codes in a straight-forward manner. Knowing 
the viscous parameters would also provide criteria for tuning Boltzmann algorithms so that they are consistent with 
quantum transport considerations. The tuning could be accomplished by changing either the sampling factor or 
by adjusting scattering algorithms. The simple manner with which the viscous parameters scale with density and 
cross-section should simplify such a procedure. 

As transport theories address the first one fm/c of a relativistic heavy ion collision, the role of non-local interactions 
becomes increasingly important. For times above 2 fm/c, it is unlikely that the non- localities play any significant role 
as the effects scale as t~^ . Since the nuclei pass one another at RHIC in less than 0.2 fm/c, whereas cross sections 
might approach a square fm, non-local effects might provide a non-negligible source of stopping as the viscous drag 
converts longitudinal collective velocity to heat. The role for such effects in the stopping phase at LHC collisions 
should be even greater. 

Finally, we mention the effect of non-local interactions one might expect in the earliest moments of the big bang. 
In that case cross sections should become perturbative and particles should be approximately massless. In this case 
cross sections should scale as a^/T^, where a is the unified coupling constant. Since T would scale as 1/t, the Navier- 
Stokes viscosity and the non-local viscosity should both scale identically with r. Since the non-local terms scale as 
a^ while the Navier-Stokes terms scale as (t~^ , the terms would differ in importance by a constant proportional to 
a^ . Thus, if the system becomes perturbative, one expects the non-local terms to be negligible compared to the 
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Navier-Stokes terms. However, the physics of high-energy hadronic coUisions is far from perturbative, and the large 
couphng constants magnify the importance of the non-local terms relative to Navier-Stokes terms. 
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FIG. 1. Spectra from GROMIT-r are displayed for four cross sections. 
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FIG. 2. Spectra for pions and protons resulting from GROMIT-r with 40 mb cross sections for four sampling factors, 
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while in the lower panel each model used a sampling factor A — 16. Due to different scattering algorithms, the models generated 
different results for A = 1, while generating identical results for larger A. 
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FIG. 4. Elliptic flow is shown as a function of transverse energy for both pions and protons for several sampling factors, 
A. The cross sections were not varied as A was varied, but instead only a fraction 1/A of the particles were scattered. The 
insensitivity to A demonstrates that the Boltzmann limit is effectively realized, even for A = 1. 
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FIG. 5. Again, elliptic flow is shown as a function of transverse energy for both pions and protons for several sampling 
factors. In this case the cross sections were scaled by a factor 1/A to account for the over-sampling. The reduced flow resulting 
for calculations with small A (larger cross sections) is attributed to viscous effects arising from the flnite interaction range. 
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FIG. 7. The mean transverse energy for particles participating in a boost-invariant Bjorken expansion is displayed as a 
function of the proper time (circles). Also displayed are the Bjorken hydrodynamic result (dotted line), the Navier-Stokes 
correction which accounts for viscous shear arising from a finite mean free path (dashed line) and the correction due to 
non-local interactions as expected from simple scaling arguments (solid line). Non-local interactions are important at small 
times when the velocity gradients and collision rates are high. 
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FIG. 8. The mean transverse energy is plotted as a function of time for the case where a forward-peaked cross section as 
described in Eq. (p3|) is implemented. The results from GROMIT-r (circles) are well described by the Navier-Stokes (dashed 
line) correction to the Bjorken solution (dotted line) . As compared to s-wave scattering, the Navier-Stokes viscosity is increased 
while the non-local contribution to the viscosity becomes negligible. 
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